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I. Introduction 


This grant to the University of Michigan supported the efforts of Michael R. 
Combi to serve as a Co-Investigator in collaboration with a larger effort by the Principal 
Investigator, William Smyth of Atmospheric and Environmental Research, Inc. The 
overall objective of this project was to analyze in a self-consistent manner unique optical 
0(lD) and NH 2 ultra-high resolution line profile data of excellent quality and other 
supporting lower-resolution spectral data for the coma of comet P/Halley by using highly 
developed and physically-based cometary coma models in order to determine and explain 
in terms of physical processes the actual dynamics and photochemical kinetics that occur 
in the coma. The justification for this work is that it provides a valuable and underlying 
physical base from which to interpret significantly different types of coma observations in 
a self-consistent manner and hence bring into agreement (or avoid) apparent 
inconsistencies that arise from non-physically based interpretations. The level of effort 
for the Michigan component amounted to less than 3 person-months over a planned 
period of three years. The period had been extended at no extra cost to four years 
because the Michigan grant and the AER contract did not have coincident time periods. 
An effort of somewhat larger scope was undertaken by the PI. 

The importance of the 0( *D) profiles is that they provide a direct trace of the water 
distribution in comets. The line profile shape is produced by the convolution of the 
outflow velocity and thermal dispersion of the parent water molecules with the 
photokinetic ejection of the oxygen atoms upon photodissociation of the parent water 
molecules. Our understanding of the NH 2 and its precursor ammonia are important for 
comet-to-comet composition variations as they relate to the cosmo-chemistry of the early 
solar nebula. Modeling of the distribution of NH 2 is necessary in order to infer the 
ammonia production rates from NH 2 observations. 

II. Accomplishments for the University of Michigan Effort 

This project represented a very well-focused effort to analyze interferometric 
spectra of comet P/Halley, obtained by Fred Roesler, Frank Scherb and student co- 
workers at the University of Wisconsin during 1985 and 1986. Since the technical details 
and scientific results are presented in great detail in the accompanying paper (Smyth et al. 
1994) only a brief summary of the work and the results will be given specifically as part 
of the main body of this report. This paper has already been submitted for publication in 
The Astrophysical Journal. 
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The portion of the work done at the University of Michigan involved the actual 
modeling analysis of the line profiles of 0( ] D) and NH 2 provided by the Wisconsin 
group and participation with the PI and the Wisconsin group in the interpretation of the 
data. A slight variation of our coma models (Combi and Smyth 1988 a&b; Combi 1989), 
which had been to analyze similar but already published data, including radio line 
profiles, was used for the analysis of the Wisconsin interferometric data. 

The principal results were as follows: 

1. The Monte Carlo Particle Trajectory Model provided excellent agreement with 
both O^D) and NH 2 line profiles from 1986 data obtained in January (r=0.78AU) and 
May (r=1.67AU). This verifies the underlying dynamics, exothermic photodissociative 
chemisty, and collisional thermalization in the coma. 

2. The somewhat wider 0( 1 D) line profile in January compared with May confirms 
the variation of the coma outflow speed with heliocentric distance. 

3. The narrow January NH 2 line profile is indicative of the outflow dynamics of 

the parent molecules in the coma because of the very short photodissociation time of the 
NH 2 parent (presumable NH 3 ) and as such is consistent with the dynamics implied by 

the O(^D) line profile and other data. 

4. The average water production rates derived from the O^D) line profiles were 
2.90 and 2.68 x 10 30 molecules per second for January 16 and 17, respectively. They 
were found to be consistent with the extrapolated water production rates determined by 
the previously published low resolution observations of O(^D) and Hot (Smyth et al. 
1993) and the 7.6-day photometric light curve of Schleicher et a. (1990). These large 
values establish that the maximum water production rate for comet Halley actually 
occurred during the 2-3 weeks prior to perihelion. Implications drawn from the 
collisionally quenched OH 18-cm radio observations in later January indicate that the 
water production rate may have been as large as 3.6 x 10^® molecules per second near 
January 23, but this is more speculative. 

5. The NH 3 production rate determined from the NH 2 line profile of January 17 
was 1.48 x 10 28 per second, yielding an NH 3 /H 2 O production rate ratio of 0.55%, which 
is consistent with most other measures from NH 2 and NH observations when analyzed 
with improved and consistent g-factors for NH 2 - 

6 . Perhaps most importantly, the MCPTM in combination with the hybrid dusty- 
gas dynamic/Monte Carlo model have now been successful in analyzing and explaining a 
wide variety of observations of the water dissociation products in comet Halley, including 
line widths of published radio OH profiles (Combi 1989), wide field Lyman-a images 
and scans (Smyth, Combi and Stewart 1991; Smyth et al. 1994), IUE Lyman-a 
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observations (Combi and Feldman 1993), IUE OH observations (Combi, Bos and Smyth 
1993), and O(^D) and Ha observations (Smyth et al. 1993). These results now present 
reasonably self-consistent analyses of water production rates and water dissociation 
product kinetics and dynamics. 
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ABSTRACT 


A set of high resolution Fabry-Perot measurements of the coma of comet 
P/Halley were acquired in the [O I] 6300 A and NH 2 6298.62 A emission lines. 
These high resolution measurements provide the first optical observations capable 
of studying directly the photochemical kinetics and dynamic outflow of the coma. 
The observations were analyzed by a Monte Carlo Particle Trajectory Model. The 
agreement of the model and observed line profiles was excellent and verified the 
underlying dynamics, exothermic photodissociative chemistry, and collisional 
thermalization in the coma. The somewhat wider line profile width for the 0(*D) 
emission in January 1986 compared to May 1986 is, for example, produced by the 
larger outflow speed nearer perihelion in January. The January NH 2 profile, 
which is more narrow than the January O^D) profile, is indicative of the outflow 
dynamics of the parent molecules in the coma. The absolute calibration of the 
observations in January 1986 allowed the production rates for the H 2 O and NH 2 
parent molecules to be determined. The average daily water production rates 
derived from the 0(*D) emission data for January 16 and 17 are, respectively, 
2.90 x 10 30 molecules s’ 1 and 2.68 x 10 30 molecules s* 1 . These very large water 
production rates are consistent with the extrapolated (and 7.6-day time variable) 
water production rates determined from the analysis of lower spectral resolution 
observations for O^D) and Ha emissions (Smyth et al. 1993) that covered the 
time period up to January 13. The large production rates on January 16 and 17 
establish that the maximum water production rate for comet Halley occurred pre- 
perihelion in January. Implications drawn from comparison with 18-cm radio 
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emission data in January suggest that the peak water production rate was even 
larger and might have been as large as ~3.6 x 10^° molecules s ' 1 near January 23. 
The average production rate for NH 3 determined from the NH 2 emission data for 
January 17 was 1.48 x 10 28 molecules s* 1 , yielding a NH 3 /H 2 O production rate 
ratio of 0.55%. This ratio is consistent with the range of earlier derived values. 
The corrected g-value noted in the NH 2 analysis brings most NH 2 production 
rates in line with NH production rates and also makes them consistent with the 
production of both species from NH 3 . 


Subject headings: Comets 



1. INTRODUCTION 


An important goal of cometary physics is understanding the 
photochemical kinetics and dynamics of molecules within the coma. Although 
radio emission lines from the coma have yielded much valuable information on 
coma gas dynamics, the possibilities of using optical emission lines for this 
purpose have not been fully explored. The relatively slow speed of cometary 
molecules leads to very narrow optical emission lines. The resolving power and 
sensitivity required for measuring the profiles of the very narrow and faint lines is 
difficult to achieve using conventional grating techniques. However, as shown in 
this paper, the problem is readily handled using Fabry-Ferot spectroscopy. 

This paper reports a unique set of high resolution Fabry-Perot 
measurements of [O I] 6300 A and NH 2 6298.62 A emissions from the coma of 
comet Halley, and its successful interpretation using advanced coma gas dynamic 
and photochemical models. It provides the first conclusive demonstration that 
high resolution optical measurements, in addition to radio measurements, are 
important and powerful means for studying the details of coma gas dynamics. 
The observations were acquired by the Wisconsin Space Physics Group in 
January and May 1986 at the McMath-Pierce solar telescope on Kitt Peak. The 
[O I] 6300 A line profiles were particularly diagnostic for studying contrasting 
cometary outflow conditions, since the Halley heliocentric distance in January 
was small (~0.78 AU) and near perihelion (0.59 AU) so that the H 2 O production 
rate was near its maximum value, while the heliocentric distance in May was 


4 



much larger (1.68 AU) so that the H 2 O production rate was smaller by 
approximately an order of magnitude. For comet Halley, these optical 
measurements are particularly important since they provide water production rates 
during a critical pre-perihelion period when these rates were near their maximum 
values, when few measurements with the exception of radio emission 
observations were made, and when the radio OH measurements are difficult to 
interpret because of uncertainties in corrections for collisional quenching at high 
water production rates. 

The purpose of the paper is to study the photochemical kinetics and 
dynamic outflow of the coma of comet Halley by analysis of the high resolution 
optical observations reported herein. A description of the instrument and 
observations is given in section 2. The cometary model and the analysis of the 
observations together with their implications are discussed in section 3. A 
summary of the results of the paper is presented in section 4. 

2. INSTRUMENTATION AND OBSERVATIONS 

The instrument used for these observations was a Wisconsin 150 mm dual- 
etalon, pressure-scanned Fabry-Perot spectrometer installed at the McMath-Pierce 
solar telescope on Kitt Peak where it could receive light either from the main 
heliostat or the west auxiliary heliostat. The instrument was used for a variety of 
observations, most of which were carried out at a medium spectral resolving 
power of 30,000 (Reynolds et al. 1986; Magee-Sauer et al. 1988, 1989, 1990; 
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Scherb et al. 1990; Smyth et al. 1993). The conversion of the spectrometer to the 
high-resolution mode (R = 190,000) for the line profile studies was effected by 
replacing one medium-resolution etalon with a high-resolution etalon. The 
general layout of the instrument and modifications are discussed in Scherb et al. 
(1986) & Roesler et al. (1986). Details important for the present paper are 
discussed below. 

For observations conducted prior to March 1986, the high resolution scans 
were taken using the main heliostat, which uses a 2.0 m flat to feed a 1.52 m 
diameter, 82.5m focal length mirror. The sky was imaged through a guider onto a 
160 mm focal length lens that imaged the primary mirror on the 3.2 mm diameter 
(slightly oversized) entrance aperture of the instrument. The light was 
recollimated for the etalon train by a 762 mm focal length lens. A field lens at the 
entrance aperture reimaged the sky through the Fabry-Perot collimating lens onto 
the midplane of the etalon train. With this coupling arrangement, the effective 
field of view was not crisply defined because of vignetting in the etalon train and 
the need to mask out a defective region of the high-resolution etalon plates. For 
the medium resolution studies, in which absolute intensity measurements were a 
primary objective, spatial scans made by letting a star drift across the field of view 
were used to map the field response and determine the effective field of view 
(Magee-Sauer 1988). This mapping was not done for the high resolution 
observations, and the map for the medium resolution mode cannot be applied to 
the high resolution mode. However the effective field of view diameter is 
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reasonably estimated to be 1.0 arcmin, and can be accurately boxed between 1.1 
and 0.9 arcmin. 

The instrumental spectral profile (ISP) was determined by scanning the 
6328 A line from a stabilized He-Ne laser, it is not expected that the profile was 
significantly different at 6300 A. This was corroborated by scans of the Cel 
6300.21 A line from a hollow cathode lamp and by thermospheric temperatures 
estimated from scans of the airglow 6300 A line in some of the comet scans. For 
all high-resolution measurements the etalon spacers were 6.25 mm and 0.472 mm, 
respectively. This choice gave appropriately high resolution, but left ghosts in the 
transmission function as high as 10% due to incomplete suppression of the 
transmission peaks of the high resolution etalon. The relative strengths of the 
ghosts were measured in a scan of the 6300 A airglow line and are included in the 
fit used to measure the cometary and airglow line widths. The observation 
periods were chosen so that ghosts of the 6300 A airglow line could be accurately 
separated from the Doppler-shifted cometary line (see Figure 1). 

After March 1986, the optical configuration was changed to accommodate 
an optional imaging detector used for some of our other programs. In the new 
arrangement, the sky was imaged on the entrance aperture of the Fabry-Perot 
system, and the telescope primary was imaged within the etalon train (see Roesler 
et al. 1986). All of the high-resolution scans with the instrument in this 
configuration were made using the west auxiliary telescope. The field was cleanly 
determined by the entrance aperture diameter, and the scale was measured both by 
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imaging a star field and by timing a star drifting a known distance across the field. 
The field of view diameter for the May 1986 high-resolution observations was 3.5 
arcmin. As before, the instrumental profile was determined with the stabilized 
He-Ne laser at 6328 A, with corroborating evidence for its applicability from Ce 
6300.21 A and airglow [O II 6300 A line widths. He-Ne laser scans were made 
with various smaller entrance aperture diameters to confirm that spatial variations 
in intensity across the field would not noticeably bias the measured profiles. 

Observational information for the high-resolution [O II 6300 A and NH 2 
6298.62 A scans is summarized in Table 1. All observations were taken with the 
field of view centered on the comet head. Selected observations are illustrated in 
Figures 1-3. In Figure 1, a high -resolution [O I] 6300 A scan obtained on 16 
January 1986 is shown. The airglow and comet (O I] 6300 A emission lines are 
the prominent features. Other apparent emission features are the combined ghosts 
from the airglow and comet. At this high resolution, there is no NH 2 
contamination of the comet [O I] 6300 A emission. Figure 2 shows a scan of NH 2 
6298.62 A emission obtained on 17 January 1986, and Figure 3 shows the sum of 
four [O I] 6300 A scans obtained on 3 May 1986. In Figure 1 and Figure 3, the 
comet emissions are clearly wider than the airglow. In Table 2, the field of view 
and the average full width half maximum (FWHM) for each day of observations 
in Table 1 are summarized. The spectral resolution of the instrument is also 
indicated in Table 2 and is more than sufficient to resolve the measured line 
profiles. The ISP used to determine the spectral resolution was obtained from 
Gaussian fits to scans of the stabilized He-Ne laser. TwoGaussians were needed 
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to fit the narrow core and broad wings of the instrumental profile. Gaussians 
convolved with the appropriate instrumental profiles were fitted to the emission 
line data in order to determine the FWHM of the measured line profiles given in 
the last column of Table 2. 

Although absolute intensity measurements were not originally planned as 
part of the high-resolution observations, a scan of the hydrogen Ha line, taken on 
the Trapezium region of the Orion Nebula on UT 18 January 1986, using the high 
resolution mode, provides a means of extracting intensities from the Halley [O I] 
6300 A and NH 2 scans taken on 16 and 17 January 1986. The basis for this 
procedure is an intensity-calibrated, CCD Ha image of the Orion Nebula obtained 
by Dufour and Hester at Palomar Observatory (Walter, Dufour & Hester 1992). 
The image is 16 arcmin square with a resolution of 1.192 arcsec per pixel (Walter 
1992). We found that the Ha intensities in this image are in good agreement with 
Ha intensity results from an imaging spectrophotometric study of the Orion 
Nebula by Pogge, Owen & Atwood (1992, 1993). Corrections for variation of 
instrumental sensitivity at the different wavelengths of Ha and the comet lines 
were obtained from earlier observations of the star a CMi (Magee-Sauer et al. 
1990; Magee-Sauer 1988). Atmospheric transmission corrections were based on 
average transmission conditions at Kitt Peak. The Orion scan data were corrected 
for the effects of the ghosts in the instrumental profile. The overall uncertainties 
in the comet intensities given in Table 1 are estimates based on uncertainties in 
the absolute Orion Nebula Ha intensity, atmospheric transmission corrections, 
and corrections for wavelength dependence of instrumental sensitivity. 
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3. ANALYSIS OF THE OBSERVATIONS 


3.1 Model Description 

In order to calculate the high resolution line profiles for both the [O I] 
6300 A and NH 2 6298.62 A emissions, we have used the Monte Carlo Particle 
Trajectory Model (MCPTM). For NH 2 and for the OH component of the [O I] 
6300 A emission, the most general heavy species collision algorithm as recently 
presented by Combi, Bos & Smyth (1993) has been used. This MCPTM is the 
most recent update to the detailed description provided in previous papers by 
Combi & Smyth (1988a, b) which had its origins in the first Monte Carlo model 
by Combi & Delsemme (1980). In the MCPTM, the description adopted for the 
outflowing collisional inner coma was determined from a hybrid gas- 
dynamic/Monte Carlo model calculation for H atoms produced photochemically 
from the H 2 O exothermic dissociative chain. In this hybrid calculation, the Monte 
Carlo model for H atoms was used to compute the photochemical heating 
efficiency of the coma and was run in an iterative fashion with the dusty-gas 
dynamic model in order to compute consistendy the gas outflow and temperature 
of the coma for the conditions expected in both the January and May time periods. 
In the hybrid calculation, a reasonable estimate of the water production rate is 
required and was obtained from many previously published data. Table 3 shows 
the model parameters used in the hybrid model calculation. These parameters and 
the water production rate used are equivalent to those adopted by Combi (1989) 
and more recently used in the analysis of the February 2-6, 1986, Pioneer Venus 
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Orbiter UVS image of the H coma of Halley (Smyth, Combi & Stewart 1991). 
For the purposes of the basic coma modeling, water production rates were used 
which are accurate to within 20-30%. Differences of a factor of two or more 
begin to have noticeable effects on the coma velocity and temperature Fields, 
whereas a factor of ten has a large effect (Bockelee-Morvan & Crovisier 1987; 
Combi 1989). The hybrid model calculations for the variations of the outflow 
speed and the gas temperature, with distance from the nucleus, are shown in 
Figures 4a and 4b for both the January 16-17 and May 3-4 epochs. Note the much 
larger outflow speeds and gas temperatures in January compared with May. 
These outflow speeds are consistent with various radio observations of OH 
(Schloerb, Claussen & Tacconi-Garman 1987; Bockelee-Morvan, Crovisier & 
Gerard 1990) and HCN (Schloerb et al. 1987). 

With the conditions of the inner collisional coma determined for the 
observational dates in January and May, the MCPTM calculations for the 
observed line profiles can be performed. Line profile information is saved in the 
MCPTM by accumulating the line of sight component of all atoms or radicals 
within the projected observational aperture in a set of velocity bins (Schloerb & 
Gerard 1985; Combi 1989) which were 0.1 km s' 1 in width in our case. Because 
of the speed of modem workstations and in order to obtain sufficiently high 
statistical accuracy, each of the models presented in this paper used 5,000,000 
simulation molecules. Full collisional models in a coma of high gas production 
take only a fraction of an hour to run. Values adopted in the MCPTM for the 
photodissociation lifetimes of H 2 Q and OH and for their corresponding branching 
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ratios that determine the O('D) production rate are the same as discussed by 
Smyth et al. (1993). For NH 2 , lifetimes were adopted from the results of the 
spatial profile analysis of Fink, Combi & DiSanti (1991), who found Haser model 
parent and daughter scale lengths which were reasonably consistent with 
production from NH 3 photodissociation. The lifetimes used in the model, reduced 
to a heliocentric distance of 1 AU, are given in Table 3. The NH 2 lifetimes are an 
'average' of the pre- and post-perihelion values found by Fink, Combi & DiSanti, 
who found a clear asymmetry in both the parent and daughter scale lengths for 

nh 2 . 

There exists the possibility, even the likelihood, that when NH 2 is 
produced in the photodissociation of NH 3 that there will be some excess energy 
which would be available for translational velocities for the fragments, NH 2 and 
H. Since we have data only for the January, when the water production rate was 
quite large, any excess translational velocity of NH 2 would be nearly immediately 
collisionally quenched. Such behavior has been illustrated by Combi, Bos & 
Smyth (1993) even for OH which has a much longer parent lifetime and for a case 
of an even smaller water production rate. We have studied the importance of 
collisional quenching using the MCPTM and test ejection speeds for NH 2 as large 
as 1 km s'*. There was no measurable effect for ejection speeds for NH 2 . This 
may not be the case at large heliocentric distances and/or lower water production 
rates. Tegler (1992) has suggested a translational ejection for NH 2 on the order of 
0.5 km s'l 
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3.2 Line Profile Analysis 

The appearance of line profiles from various types of simple coma outflow 
models has been recently discussed in a paper by Bockel6e-Morvan, Crovisier & 
Gerard (1990). They show, for example, that a classical ’vectorial' type model 
produces a trapezoidal shaped line profile. Combi (1989) showed line profiles 
resulting from more physically realistic modeled comae that included the intrinsic 
variation of outflow speed with increasing distance from the nucleus and the gas 
temperature. These effects tend to make the unusual looking profiles appear more 
Gaussian. Real observed profiles are furthermore convolved with the ISP, that is 
here approximated by two Gaussians, and that under the best of circumstances has 
a width nearly as large as the intrinsic width of the coma line in question. For the 
analysis presented here, the ISP was determined from scans of the helium-neon 
spectral laser and was convolve with coma model line profiles. Nonetheless, it is 
instructive to examine the shape of model line profiles before convolution with 
the ISP. For the 0( l D) profile it is also useful to examine the appearance of the 
H 2 O and OH sources of 0( ! D) atoms separately. 

Figure 5a shows the O(’D) emission line profile corresponding to the 
average conditions for January 16-17 before convolution with the ISP. Figure 5b 
shows a similar plot for May 3-4. Shown are the separated contributions from the 
H 2 O and OH sources as well as the sum, weighted according to the branching 
ratios given in Table 3. In both cases the H 2 O source dominates the emission 
within the aperture, because in both cases most of the water in the coma is within 
the aperture but most of the OH is outside. The aperture was larger in May 
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compared with January, but proportionally smaller compared with the much larger 
photodecay lifetime of water, so that the relative OH-source contribution was 
even smaller, as shown. If the entire coma were observed, the OH source would 
actually be nearly twice that of the H 2 O source. The profile from the OH source 
is somewhat wider than that from H 2 O because of the additional 1.05 km s * the 
OH receives upon dissociation from water. This is counteracted slightly by the 
ejection speed of the O^D) atoms which is actually slightly smaller from OH than 
from H 2 O (1.5 vs. 1.6 km s'l). 

Figures 6 a&b show an example of the effects of the ISP convolution on 
the modeled profiles for O(^D) and NH 2 in January. The slight apparent offset of 
the intrinsic cometary and ISP convolved profile is due to the fact that the actual 
ISP from the He-Ne laser was slightly asymmetric. The ISP causes the final 
observed profile to appear very Gaussian like, and removes any peculiar or 
interesting features that the profile could show if observed at a substantially 
higher spectral resolution. 

The agreement between modeled and observed line profiles is excellent. 
Figures 7 a&b show observed 0( ! D) profiles from January and May compared 
with the model profiles for the combined H 2 O and OH sources convolved with 
the ISP. For the May data we have summed all of the profiles observed on May 3 
and 4 to improve signal-to-noise. The somewhat wider January profile is 
naturally explained by the larger outflow speed as predicted by the gas- 
dynamic/Monte Carlo model (Figures 4 a&b). Figure 8 shows the observed NH 2 
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profile from January 17 compared with the comparable model again convolved 
with the ISP. Because of the combination of the high water production rate and 
short lifetime for the NH 2 parent (NH 3 ), the width and shape of the NH 2 profile 
should be indicative of the outflow speed and temperature for parent molecules 
(H 2 O). As mentioned above, especially for the January coma conditions, we 
expect there to be no measurable effect from any photodissociative ejection 
velocity for NH 2 , if one exists, because of collisional quenching. The slight 
misalignment (< 0.2 km s ’ 1 or < 4 mA) of the center of the model and observed 
line profiles in Figures 7 a and 8 (i.e., the model profiles need to be moved slightly 
to the left in the January data) and in Figure 7b (i.e., the model profile needs to be 
moved slightly to the right in the May data) is likely due to small motions typical 
in the earth's atmosphere. These motions cause the center of the terrestrial 0( ! D) 
line profile, used to determine the absolute wavelength scale for the observations, 
to vary slightly. 

In order to investigate the sensitivity of the profiles to changes in model 
parameters, we have also compared a number of simple Monte Carlo/vectorial 
models, i.e. models with single values of the parent molecule and ejection speeds 
for the daughters, for a range of outflow speeds, and convolved them with the 
ISP. It is clear that outflow speed differences more than 0.2 km s " 1 are quite 
significant and that such differences between the model and measurement would 
be easily detectable in data-model comparisons. The outflow speed differences 
between January and May are larger than this (see Figures 4 a&b). 
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3.3 Production Rates 


In addition to the line profile information, as mentioned above, an absolute 
calibration was possible for all January observations in Table 1 . With this 
information, it was possible to use the model results to calculate NH3 and H2O 
production rates. The MCPTM simulates the coma by using a large number of 
molecules (5 million in this case for good statistical coverage) that are ejected into 
the coma over a long enough time (2 x 10^ seconds) so that the model coma 
relaxes to a realistic steady-state condition. All of the computed quantities in the 
model (e.g., abundances) are normalized to the simulation production rate, which 
is the number of simulation molecules divided by the total simulation run time. 
Simulation quantities are calculated per unit production rate and can then be 
converted to physical quantities by multiplying by the actual production rate. For 
any decaying species, produced at a constant rate, the sum of all particles in the 
coma is simply the product of the production rate and the lifetime of the species. 
Since the model is normalized to a unit production rate, the total number of 
simulated molecules (or radicals) in the coma is equal to the lifetime of that 
species. In the models for the analyses presented here, the total number of 
simulation radicals within the aperture is obtained simply by integrating over the 
line profile. When this number is multiplied by the production rate, the number of 
actual radicals within the aperture is obtained. 

For the case of NH2, the total number of radicals in the aperture, times the 
fluorescence rate (or g-factor), in photons per second per radical into 4 ji 
steradians, times 10'^, yields the intensity in Rayleighs. Since the brightness is 
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observed, the production rate is the only unknown and can be calculated. The 
absolute value of the fluorescence rates for the various bands of NH 2 has been a 
contentious issue in the last several years. In the pre-Halley review paper by 
A'Heam (1982), the g-factors for the (0,7,0) through (0,12,0) NH 2 bands were 
summarized as computed in an earlier paper (A'Heam, Hanich & Thurber 1980). 
This was the result of a simple calculation assuming that fluorescence occurs in 
each band individually from the ground state. Tegler & Wyckoff (1989) 
performed a more sophisticated calculation extending to more bands, (0,3,0) 
through (0,15,0), and using updated experimental data and theoretical transition 
probabilities. The newer g-values were smaller by factors of 1.4 to 5.9. Most 
importantly for our purposes here, the g-value for the ( 0 , 8 , 0 ) band decreased from 
1.694 x 10 * 2 to 4.95 x 10 * 3 photons per second. Such a change would increase 
NH 2 production rates determined from (0,8,0) band observations by a factor of 
more than three. Finally, the "Note added in proof' in the paper by Magee-Sauer 
et al. (1989), quotes Arpigny as finding that there are Ka"=l levels which are not 
sampled by the (0,8,0) band. This results in "hidden" NH 2 , if one were to 
determine NH 2 abundances from (0,8,0) band observations alone. The effect is 
about a factor of 2 , which would bring the g-value down even further to 
2.475xl0 * 3 photons per radical per second into 4 k steradians. This extra factor 
of 2 brings most NH 2 production rates in line with NH production rates 
(Schleicher & Millis 1989; Feldman et al. 1993) and makes them consistent with 
production of both species from NH 3 . Using this g-value, the NH 2 production 
rates determined from the 17 January 1986 high-resolution observations are 
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summarized in Table 4. The average NH 2 production rate for 17 January was 
1.48 x 1028 molecules s'*. 


For the case of 0( ! D), the calculation of a water production rate is more 
complicated, but still relatively straightforward. Since the O^D) emission is a 
spin-forbidden transition it does not emit by fluorescence but rather as a result of 
photodissociation of a fraction of H 2 O, producing H 2 +0( ! D), and a fraction of 
OH, producing H + 0( ! D). The best estimate for the branching ratio from the 
H 2 0 source is from Smyth et al. (1993) for the quiet sun conditions relevant for 
the Halley time period. This value is consistent with the earlier estimate of Combi 
& Smyth (1988b). Since the OH lifetime is dominated by predissociation in the 
near UV (Schleicher & A’Heam 1988) and is highly dependent on the heliocentric 
velocity of the comet, the branching ratio actually varies even though the rate for 
the far UV component that is responsible for the 0( l D) production (van Dishoeck 
& Dalgamo 1984) does not vary with heliocentric velocity. For January 16 and 
17 the average heliocentric velocity was -23.3 km/s. The values of all of the 
relevant branching ratios and lifetimes are given in Table 3. 

The brightness of 0( ! D) per unit production rate was calculated from the 
MCPTM given the total abundances of H 2 0 and OH within the aperture (which 
were obtained by integrating over the line profiles) times the total production rate 
of 0(*D) atoms from each parent species times 10‘ 6 to get brightness units of 
Rayleighs. These rates, which are defined in Table 3, are given as 
v(H 2 0)=(BR2)/ti and v(OH)=(BR1)(BR3)/x 2 . Finally, one needs to account for 
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the 3 to 1 ratio between the two lines of the *D to 3 P doublet at 6300 and 6364 A, 
respectively. This introduces a factor of 3/4 of a 6300A photon emitted per O('D) 
atom produced. The calculated H 2 O production rates determined in this manner 
from the 0( J D) observations arc given in Table 5. The water production rates in 
Table 5 are quite high. The average value on January 16 is 2.90 x 10 30 molecules 
s"l, and the single value on January 17 is 2.68 x 10 3 ® molecules s - *. These values 
are in fact the highest water production rates yet reported for comet Halley, higher 
even than those right around perihelion (Smyth, Combi & Stewart 1991; Smyth, 
Marconi & Combi 1993). 

3.4 Halley's Activity in January 1986 

Because of its close proximity to the sun in the sky during January 1986, 
most ground-based optical observations of the comet stopped by January 12 or 13 
(Schleicher et al. 1990, Fink, Combi & DiSanti 1991, Smyth et al. 1993). 
Therefore, there are little directly comparable data with which to compare our 
production rates from January 16 and 17. In order to provide a check on the 
validity of these large water production rates, we have, however, compared these 
values with the published medium spectral resolution values from January 4 
through 13 recently reanalyzed by Smyth et al. (1993) and shown to be consistent 
with a comparable set of H Balmer-a observations. In addition, we have also 
compared those early January water production rates with a nearly simultaneous 
set of C 2 production rates to examine the correlation of day-to-day changes in 
production rates between species and how the expected 7-plus day periodicity in 
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those data (Schleicher et al. 1990) might extrapolate to the production rates on 
January 16 and 17. 

To facilitate this extrapolation, it is important first to compare the earlier 
water and C 2 production rates. Figure 9 shows the production rates of water from 
January 4 through 13 in the analysis of medium-spectral-resolution 0( ! D) and Ha 
observations of Smyth et al. along with the C 2 production rates during the same 
period published by Schleicher et al. and with new unpublished C 2 rates very 
recently obtained from Schleicher (1993) for January 11, 14 and 18 (i.e., days 
from perihelion of about -29, -26 and -22). These three points were corrected so 
as to be consistent with the C 2 scale lengths used previously by Schleicher et al. 
(1990). From day -38 to -27 where the data coverage is almost daily, it is clear 
that the same general day-to-day time variation of the C 2 production rate is 
evident in the water production rate determined from the 0( ! D) and Ha 
observations. The photometry was not abundant enough for Schleicher et al. to 
derive a very good period for January 1986 alone. When grouping all C 2 data 
from November 1985 through January 1986, Schleicher et al. found an optimal 
value of 7.46 days. This is in contrast to 7.60 days for their March data (which 
included late February through mid-March ) and 7.37 days for their April data 
(which included late March through late April). Combi & Fink (1993), 
furthermore, showed that the analysis of spatial profiles from March 1 and 2 — 
during a data gap in the photometry — required the use of the 7.60 day period as 
contrasted with the 7.37-day period. They also noted that the day-to-day 
variations of the water production rates as calculated from the Pioneer Venus 
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Orbiter observations (Stewart 1987) clearly showed that the same 7.6-day 
"March-lightcurve" was repeating at least as far back as the time of perihelion 
(February 9). 

The daily averaged water production rates for January 16 and 17 in Figure 
9 (filled circles) lie only slightly higher than the trend set by the early January 
water production rates noted in Figure 9. If the C 2 production rates are multiplied 
by 832, they match the water production rates very well from day -38 to -27 
where the data coverage is almost daily. An extrapolation of these rates for later 
times was then accomplished by shifting them ahead by 7.60 days and adjusting 
them upwards by the heliocentric distance trend (r 3 - 25 ) of the C 2 data. This 
extrapolation is shown in Figure 10 where the water production rates determined 
from the Ha data have not been included for clarity. The January 16 and 17 
points are seen to be only slighdy higher than the extrapolated pattern and have 
the proper relative values with respect to the time dependent 7.6 day extrapolated 
pattern seen a week earlier around January 8 and 9. The high resolution point on 
January 16 falls just after the peak seen in the extrapolated medium resolution 
data on January 8. Examination of the individual C 2 light curves for different 7- 
plus day periods (Schleicher et al. 1990) shows two different time patterns for the 
maximum amplitude, the larger of which is very similar to the pattern produced 
by the two water production rates on January 16 and 17. We conclude that the 
water production rates on January 16 and 17, although large, are very reasonable 
and have values similar to those expected from extrapolating earlier data. 
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A comparison of the average H 2 O production rates derived from the high 
resolution O^D) observations on January 16 and 17 with other H 2 O production 
rates is presented in Figure 11. As noted in Figure 9, the water production rates 
on January 16 and 17 (• symbols) appear to be a continuation of the rates 
determined from the O^D) data (Q symbols) and Ha data (V symbols) of 
Smyth et al. (1993) which were acquired up to January 13. From this comparison, 
the importance of the pre-perihelion peak in the water production rate can be 
appreciated even though there is a gap in the daily coverage (with the exception of 
the radio emission data, x symbols) from January 18 to January 31, where the 
Pioneer Venus H Lyman-a measurements (+ symbols) then begin on February 1. 
The largest production rate near perihelion determined from the Pioneer Venus H 
Lyman-a measurements is 1.8 x 10 3 ^ molecules s**, which is a factor of 1.6 
smaller than the 2.90 x 10 30 molecules s' 1 value determined for January 16. The 
water production rate was, however, likely even larger in later January since its 
rapid rise, which in Figure 1 1 began about January 10 and is still evident on our 
last data point on January 17, shows no change in its steep upward slope. The 
water production rates derived from the 18-cm radio emission data (x symbols) 
fill in this gap, but are much smaller than the water production rates determined 
from optical and ultraviolet observations on either side of the gap by about a 
factor of 3. The factor of 3 is most likely caused by uncertainties in correcting the 
radio data for collisional quenching at these very high water production rates. If 
the water production rates for the radio emission data in the region of the gap are 
multiplied by a factor of three, their steady increase and then decrease in this gap 
would imply that the water production rate reached a pre-perihelion maximum 
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near January 23 of about 3.6 x 10 30 molecules s' 1 . Further efforts to improve the 
collisional quenching correction of the radio emission data are clearly important 
to refine this pre-perihelion behavior for the water production rate. In additional 
to these O^D) results, our reasonable values for the implied NH3/H2O production 
rate ratio (see below) places more certainty on the surprisingly large value of the 
pre-perihelion water production rate. 

3.5 The NH3/H2O Production Rate Ratio 

The yield of NH2 from NH3 photodissociation is about 95%. Most of the 
NH2 is subsequently photodissociated and yields mostly NH, which is also 
observed in comets. There has been some controversy in the last couple of years 
regarding NH3 production in comets in general and in Halley, in particular. The 
major discrepancies have been due mostly to uncertainties and overestimates of 
the NH2 g-factor, although there is some degree of uncertainty from the NH 
lifetime and/or parent and daughter Haser scale lengths which are used in the 
analysis of NH observations. Recently, Arpigny et al. (1993) reported that most 
estimates of the NH3/H2O ratio from both NH and NH2 observations in Halley 
are consistent with values somewhat less than but near 1%. Given the latest and 
best estimates for both the CK'D) branching ratios and the NH2 g-factors, our 
value for the NH3/H2O ratio on 17 January is 0.55%. The average value 
determined from the medium resolution data of Magee-Sauer ct al. (1989) for 
January 1986, corrected for the Tegler & Wyckoff (1989) band emission rates and 
the factor of 2 from Arpigny (see note in Magee-Sauer et al.), was 0.75%. 
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Other determinations for the NH3/H2O ratio from NH2 observations in 
Halley are 0. 1-0.7% by Wyckoff et al. (1988) and 0.3% by Fink (1993). All these 
have been renormalized to the same g-factor. A value determined from NH 
observations is 0.44-0.94% by Feldman et al. (1993). Schleicher & Millis (1989) 
have published a value of 0.5% for a 'normal comet as determined from aperture 
photometry of comets P/Halley, P/Borrelly & Liller (1988a). Feldman et al. 
(1993) present values of the NH/OH ratio in a number of comets observed with 
IUE to be in the range of 0.39-0.95%. These would imply NH3/H2O ratios in the 
range of 0.36-0.88%. Finally, all of these 'measured' NH3/H2O ratios depend on 
the models and model parameters employed in their respective analyses; the 
analyses are far from uniform. Despite the nonuniformity of analysis, it is likely 
that both NH and NH2 observations are consistent with production from NH3 in 
the range of 0.4 to 1.0% of that of water. Our value of 0.55% fits well within this 
range. 

4. SUMMARY 


The first high spectral resolution optical observations for the coma of 
comet Halley capable of studying the photochemical kinetics and dynamics of 
outflowing molecules were acquired in 1986 for [O I] 6300 A and NH2 6298.62 A 
emissions in the January pre-perihelion time frame and also for [O I] 6300 A 
emissions in the May post-perihelion time frame. The measured line profiles 
were obtained for apertures centered on the nucleus, and an absolute calibration of 
the emission line brightnesses in January was determined. The observations were 
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analyzed by employing the MCPTM of Combi & Smyth (1988a, b), using the 
recent heavy species collisional algorithm of Combi, Bos & Smyth (1993) for 
NH2 and the OH component of [O I] 6300 A emission. 

The agreement of the model line profiles with the observations in January 
of 0( ! D) and NH2 and in May of 0( 1 D) was excellent and represents yet another 
in a series of verifications in the last several years of die modeling techniques and 
physical parameters employed. Model simulations show that outflow speed 
differences more than ~0.2 km s" 1 are quite significant and that such differences 
between the model and measurement would be easily detectable in data-model 
comparisons. The somewhat wider line profile width for the O(^D) emission in 
January compared to May is produced by the larger outflow speed nearer 
perihelion in January. Doppler broadening of the NH2 line profiles because of the 
(expected, but yet undetermined) translational velocity imparted to NH2 fr°m its 
exothermic NH3 photodissociation was shown in model simulations to be 
collisionally quenched for molecules in the observing aperture (i.e., 0.5 arcmin 
from the nucleus) and hence unimportant. For the O(^D) line profile, the 
translational velocities imparted to 0(*D) upon exothermic dissociation of H2O 
and OH are not collisionally quenched and do play a role in the line profile width. 
The separate contributions to the O('D) emission line profile produced by the 
photodissociation of H2O and OH were calculated and discussed, and the 
composite line profile was constructed by their properly weighted sum. The 
profile contribution from OH is somewhat wider than from H2O because of the 
additional 1.05 km s' 1 that OH receives upon dissociation from water. The 
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convolution of the composite modeled line profile with the instrumental spectral 
profile, however, causes the final observed profile to appear very Gaussian like 
and removes any interesting features that the profile could show if observed at a 
substantially higher spectral resolution. 

The absolute calibration for the observations in January allows the 
production rates for the H2O and NH3 parent molecules to be calculated. The 
average water production rates derived from the O('D) emission data for January 
16 and 17 are, respectively 2.90 x 10 30 molecules s l and 2.68 x 10 30 molecules s' 
l 4 These very large water production rates are consistent with the extrapolated 
(and 7.6-day time variable) water production rates determined from the analysis of 
medium spectral resolution observations for 0(*D) and Ha emissions (Smyth et 
al. 1993) that covered the time period from January 4-13, 1986. Together with the 
water production rates derived from the analysis of medium spectral resolution 
observations, these large production rates on January 16 and 17 establish that the 
maximum water production rate for comet Halley occurred pre-perihelion in 
January. The previously largest water production rate determined from the H 
Lyman-a observations for times very near perihelion was 1.8 x lO 3 ^ molecules s 
1, which is a factor of 1.6 times smaller than the value on January 16. Although 
there are difficulties (as illustrated in this paper) in determining the absolute water 
production rate from the 18-cm radio emission data acquired in January and early 
February because of correcting for collisional quenching, the peak water 
production rate may be ~3.6 x 10 3 ® molecules s‘^ and occur near January 23 if 
they are approximately scaled by a factor of three using optical and ultraviolet 
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data as a guide. However, to do this scaling properly, a careful reassessment of 
the collisional quenching calculation would need to be done for the cometary 
plasma conditions appropriate for small heliocentric distances and large 
production rates. The average production rate for NH3 determined from the NH2 
emission data for January 17 was 1.48 x lO 2 ^ molecules s' 1 and produced a 
NH3/H2O production rate ratio of 0.55%. This ratio is consistent with the range 
of earlier derived values and therefore provides independent support for the high 
water production rates derived from the O('D) emission data in January. The 
corrected g-value noted in the NH2 analysis brings most NH2 production rates in 
line with NH production rates and also makes them consistent with the production 
of both species from NH3. 
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Log of High-Resolution Comet Halley Scans 
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Table 2 

Observational Information for the Wisconsin High Resolution 0(*D) and NH 2 Line Profiles 
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Table 3 

Parameters for the Hybrid Gas-dynamic/Monte Carlo Coma Model 

January 16 17 Gas Production Rate = 1 x 10^0 s"* 

May 3-4 Gas Production Rate = 1 x 10^9 s*l 
Radical-H20 collisional cross section = 3.5 x 10"^ cm '2 
H2O lifetime: xj=8.2 x 10^ seconds at 1 AU 

OH lifetime: T3=1.23 x 10 5 seconds at 0.78 AU (January 16-17 average) 

13=4.07 x 105 seconds at 1.68 AU (May 3-4 average) 

NH3 lifetime: 7.0 x 10^ seconds at 1AU 
NH2 lifetime: 5.0 x 10^ seconds at 1 AU 

Branching Ratios: BRl=OH from H2O = 0.88 

BR2=0( 1 D) from H2O = 0.034 
BR3=0( 1 D) from OH = 0.081 
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Table 4 

Production Rate of NH 2 from 6298.6 A Line Profiles 


UT Date of 
Observation 

Start 

Time 

(UT) 

Average 

Intensity 

(Rayleighs) 

Production 
Rate of NH2 
(10^8 molecules s'l) 

1986 January 17 

0204 

990 ± 140 

1.6210.23 


0215 

940 ± 130 

1.53 + 0.21 


0219 

800± 110 

1.29 + 0.18 


0223 

900 ±130 

1.4710.21 
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Table 5 

Production Rate of Water from [O I] 6300A Line Profiles 


UT Date of 
Observation 

Start 

Time 

(UT) 

Average 

Intensity 

(Rayleighs) 

Production 
Rate of Water 
(1030 molecules s'*) 

1986 January 16 

0132 

2240 1310 

2.8610.40 


0141 

22701320 

2.91 10.41 


0150 

22601320 

2.9010.41 


0200 

23201320 

2.9510.41 


0210 

22101310 

2.8210.39 


0215 

2200 1 310 

2.8110.39 


0221 

2200 1 310 

2.8110.39 


0225 

2430 1 340 

3.1010.43 


0230 

2390 1 330 

3.0710.43 


0235 

2170 1 300 

2.7710.39 

1986 January 17 

0228 

21001290 

2.68 1 0.38 
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Figure Captions 


FIG. 1. Scan of [0 I] 6300 A Emission, 1986 January 16, 0132 UT. The 
cometary and airglow emission line profiles are well separated and spectrally 
resolved. The airglow ghosts fall on either side of the comet [O I] 6300 A 
emission. The location of the zero on the velocity scale is arbitrary and has been 
chosen to be at the center of the cometary emission line. 


FIG. 2. Scan of NH2 6298.62 A Emission, 1986 January 17, 0204 UT. The 
cometary NH 2 emission line profile is spectrally resolved and centered on the zero 
of the velocity scale. 

FIG. 3. Sum of Four Scans of [O I] 6300 A Emission, 1986 May 3, 0204-0226 
UT. The cometary and airglow emission line profiles are well separated and 
spectrally resolved. Relative to the airglow emission line, the cometary emission 
line is much weaker in May than in January 1986 because of the smaller H 2 O 
production rate of the comet and the larger field of view of the instrument. 

Fig. 4. Results of the Gas-dynamic/Monte Carlo Model for Comet P/Halley 
in January and May 1986. Shown are the variations with distance from the 
nucleus of the coma outflow speed in (a) and of the gas temperature in (b) for 
conditions in mid-January and early May 1986. Curves are cut off at the 
reasonable limits of the hydrodynamic calculation. Note the much higher 
velocities and gas temperatures in January owing to both the higher heating rate 
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and the higher degree of collisional thermalization of the hot superthermal H 
atoms that are largely responsible for the photochemical heating of the coma. The 
effects of the higher velocities and temperatures are clearly seen even in the 
Gaussian widths of the O('D) lines given in Table 2. 

Fig. 5. Model for the H 2 O and OH Source Components of the 0( 1 D) Line 
Profile. Model results for the intrinsic O(^D) cometary line profiles for 
conditions in January are shown in (a) and those in May are shown in (b). In both 
cases the lowest profile is from the OH source, the middle is from the H 2 O source 
and the highest is the appropriately weighted sum. The instrument spectral profile 
(ISP) has not been convolved. Note that the OH source is somewhat wider than 
the H 2 O source, but that the H 2 O source is dominant for these nucleus-centered 
apertures. The line profiles were binned in 0.1 km s' * intervals. 

Fig. 6. Convolution of Intrinsic Cometary Line Profiles with ISP. The effect 
of the convolution of the ISP with the intrinsic cometary line profile shapes is 
illustrated. Shown in (a) is the Of^D) line and in (b) the NH 2 6298.6 line for 
January 1986 conditions. In both cases the lines with the higher peaks, 
irregularities and more "interesting" shapes, are the intrinsic cometary line 
profiles. The ISP tends to broaden a line and make it appear more Gaussian. In 
each case the lines have been normalized to the same integrated line flux. A slight 
asymmetry between the lines is caused by the slightly asymmetric ISP. 
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Fig. 7. Model-Data Comparison for the [O I] 6300 A Line Profile. Shown in 
(a) is the observed line profile for January 17 (the histogram) plotted with the ISP- 
convolved model line profile (smooth curve) and normalized to one another 
according to the integrated line flux. The feature centered at about +5 km s'l is a 
ghost in the data which has been properly removed in model-data normalization. 
Shown in (b) is a similar comparison for the average of the May 3 and 4 
observations. See text for a detailed discussion. 

Fig. 8. Model-Data Comparison for the NH 2 6298.6 A Line Profile. The 
average of the four observed line profiles acquired on January 17 (the histogram) 
is plotted together with the ISP-convolved model line profile (smooth curve) and 
normalized to one another according to the integrated line flux. 

Fig. 9. Coma Activity of Comet P/Halley in January 1986. Shown are our two 
new water production rates determined from the high resolution O^D) line 
profiles (plotted as filled circles) compared with those determined from the 
medium resolution 0(*D) profiles (open circles) and hydrogen Balmer-a 
(triangles) from Smyth et al. (1993). Also shown for comparison are the C 2 
production rates x 100 (open squares) from Schleicher et al. (1990) and newer 
values x 100 (filled squares) from Schleicher (1993). It is clear that the day-to- 
day structure of coma activity was the same in water, from 0(*D), as it was in C 2 . 
There are unfortunately no corroborating data in the January 16-17 time period. 
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Fig. 10. Comparison of (he Water Production Rates from High Resolution 
Line Profiles with an Extrapolated Periodic Light Curve. The two new water 
production rates determined from the high resolution 0( ] D) line profiles (plotted 
as filled circles) are compared with the trend set by the earlier medium resolution 
values (Smyth et al. 1993) which are shifted by 7.60 days and extrapolated 
according to the heliocentric distance power law. The original 0( I D) medium 
resolution data from Smyth et al. (1993) are plotted as open circles. The original 
C 2 data are plotted by open squares and filled squares. The extrapolated 0(*D) 
data are plotted as shaded open circles, and the extrapolated C 2 data are plotted by 
shaded open squares and shaded filled squares. The two new values coincide 
reasonably well with the January 8 peak, shifted and extrapolated ahead by 7.6 
days. 

Fig. 11. Comparison of the Water Production Rates for Comet Halley. The 
production rates in units of 10^ molecules/s are shown as a function of the day 
of the year from November 1985 to May 1986 for a variety of measurements. 
The day of perihelion is shown by the vertical dotted line. The symbols refer to 
the following: ( A ) H 2 O production rates determined from IUE data for OH 
emission by Combi, Bos & Smyth (1993); (CD) H 2 O production rates 
determined from the revised analysis (Smyth et al. 1993) of the 6300A emission 
data for atomic oxygen originally published by Magee-Sauer et al. (1990); (V) 
H 2 O production rates determined from the analysis of Ha emission data by 
Smyth et al. (1993); (x) H 2 O production rates determined from 18-cm OH radio 
observations by Bockel6e-Morvan, Crovisier, & Gerard (1990); (0) H 2 O 
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production rates determined from rocket observations of H Lyman-a emission by 
McCoy et al. (1992); (*) H 2 O production rates determined from [O I] 6300 A 
emission observations of atomic oxygen by Fink & DiSanti (1990); (+) 

production rates determined from the H Lyman-ct observations of Smyth, 
Marconi & Combi (1993); (•) average daily H 2 O production rates determined in 
this paper from the [O I] 6300 A emission observations in Table 1. 
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